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ABSTRACT 

We compute models of the transmission spectra of planets HD 209458b, HD 189733b, and generic hot Jupiters. 
We examine the effects of temperature, surface gravity, and metallicity for the generic planets as a guide to under- 
standing transmission spectra in general. We find that carbon dioxide absorption at 4.4 and 15 /im is prominent at 
high metallicity, and is a clear metallicity indicator. For HD 209458b and HD 189733b, we compute spectra for 
both one-dimensional and three-dimensional model atmospheres and examine the differences between them. The 
differences are usually small, but can be large if atmospheric temperatures are near important chemical abundance 
boundaries. The calculations for the 3D atmospheres, and their comparison with data, serve as constraints on 
these dynamical models that complement the secondary eclipse and light curve data sets. For HD 209458b, even 
if TiO and VO gases are abundant on the day side, their abundances can be considerably reduced on the cooler 
planetary limb. However, given the predicted limb temperatures and TiO abundances, the model's optical opacity 
is too high. For HD 189733b we find a good match with some infrared data sets and constrain the altitude of a 
postulated haze layer. For this planet, substantial differences can exist between the transmission spectra of the 
leading and trailing hemispheres, which is an excellent probe of carbon chemistry. In thermochemical equilib- 
rium, the cooler leading hemisphere is methane-dominated, and the hotter trailing hemisphere is CO-dominated, 
but these differences may be eliminated by non-equilibrium chemistry due to vertical mixing. It may be possible 
to constrain the carbon chemistry of this planet, and its spatial variation, with JWST. 

Subject headings: planetary systems; radiative transfer; stars: HD 209458, HD 189733 



1. INTRODUCTION 

Soon after the detection of the transits of plan et HD 
209458b (ICharbonneauetalJl2000t iHenrv et all feOOO). atmo- 
sphere modelers realized that transits could provide a probe of 
the atmospheric composition of transiting planet atmosphere s 
dSeaeer & Sasselovll2000t lBrownll200U iHubbardet al] 120011) . 
Stellar light, which passes through the thin upper atmosphere 
of a transiting planet, travels through a long path length in the 
planet's atmosphere, which allows atomic and molecular ab- 
sorption features in the atmosphere to be imprinted on top of 
the stellar spectrum. This predictio n was stunningly confirmed 
by observations of HD 209458b by ICharbonneau et al l d2002) 
with the STIS instrument aboard the Hubble Space Telescope. 

This first detection of an exoplanet's atmosphere began our 
current era of exoplanet characterization, which has seen our 
knowledge of these planets expand due to Hubble, Spitzer, 
MOST, CoRoT, Kepler, and ground-based efforts. The use of 
Hubble's STIS spectrograph as the premier instrument of trans- 
mission spectroscopy was cut short by the instrument's fail- 
ure in 2004, soon after the detection of an evaporating exo- 
sphere around HD 209458b (IVidal-Madiar et all 12003. 2004) 
and the full-optical t rans mission spectrum of H D 209458b 
dKnutson et al.l2007bl). ThelKnutson et all d2007bl) data set was 
further analyzed by iBarmanf (2007) a nd D. Sing and collabo- 
rators dSing et al.ll2008allbt iLecavelier des Etangs et al.l l2~008b; 
iDesert et al.ll200 8). 

For HD 209458b, the optical transmission spectra data 
sets achieved so far has basically confirmed t he predictions 
of the first m odels dSe ager & Sasselovl 120001: iBrownl 120011: 
Hub bard et al.l 12001 1) : absorption features, which can be 
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thought of as an increase in the cross-sectional area or ra- 
dius of the planet at a given wavelength, can be seen dur- 
ing the transit. Prominent features include Rayleigh scattering 
off of H2 molecules in the blue, absorption due to the neutral 
atom ic Na doublet at 589 nm, and perhaps TiO/VO absorp- 
tion dDesert et al.l2008l) and H^O absorption in the red (Barman 
2007). Zahnle et al.l d200 9b) have also pointed out that species 
such as HS may also contribute in the blue, along with the 
Rayleigh scattering contribution. 

In the absence of this optical space spectrograph, other in- 
struments aboard Hubble and Spitzer were used to probe the at- 
mos phere of planet HP 189733b. Hubble's NICMOS was used 
by I Swain et all d2008l) to obtain a near IR transmission spec- 
trum, which detected_CJi4 (but_see also [Sing^etal]|2009). Us- 
ing Hubble ACS lPont etall d2008l) detected a featureless, slop- 
ing optical spe ctrum, perhaps impl ying a Rayleigh-scattering 
haze layer, and lTinetti et al.l d2007bl) and lBeaulieu et al.l (12008) 
used Spitzer's IRAC to find absorption variation in the mid IR 
that is c onsistent with abs orption via water vapor (however, 
see also Des ert etaill2009l for a different analysis). For HD 
20945 8b at mid-infrared wav elengths, there is a 24 fim MIPS 
point (Richardson et al. 2006) and a four-band IRAC data set 
(Bea urieuetaljl2009l). HD 209458 b, a pM class planet in the 
terminology o flFortnev etaH d2008l) and HD 189733b, pL class, 
remain by far the best characterized transiting planets, and this 
will likely remain the case for quite some time. 

The characterization of these atmospheres has increased in 
sophistication, in particular with the complementary detection 
of secondary eclipses of the transiting planets, which yield 
the day-side planet-to-star flux ratios for these systems (e.g. 
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Knut son et all 120081: ICharbonneau et all 120081: iMachalek etafl 
20081) . These kinds of observations are complementary to the 
transmission spectrum observations. Due to the longer path 
length at slant viewing geometry, transmission spect ra probe 
lower atmospheric pressures at a given wavelength dFortnevI 
12005b . Perhaps more importantly, transmission spectra probe 
the terminator region of the planets, the great circle route divid- 
ing day from night. The terminator region is obviously much 
less irradiated than the substellar point or the day-side mean, 
so we could expect temperature dependent and irradiation de- 
pendent chemical mixing ratios that are different on the termi- 
nator, compared to the substellar point, due to th e lower inci- 
dent fl uxes and cooler tempera tures. Interestingly. IS wain et al.l 
(2008]) and I Swain et al.l d2009l) apparently find such a situation 
for HD 189733b, where on the day side, they find little to no 
carbon in CH4, while it is detectable at the terminator. 

More complex models of the atmospheric temperature struc- 
ture, chemical mixing ratios, and the dynamics of these at- 
mospheres have become necessary to understand these grow- 
ing data sets. Recent one-dimensional models of hot Jupiters, 
that ai med to predict the mean temperature struc t ure and spec- 
trum dBurrows et alj 120051 : iFortnev et al.l 120051 : ISeager et al.l 
l200l . and ignored dynamics, have yielded in some circum- 
stances to multi-dim ensional purely radiative-transfer models 
dBarman et al.ll2005|). or models that also include para meteriza- 
tions of dynamics dlro et al]|2005b [Burrows et al.ll2008b . 

Our group has taken a different path, first in examining 
the s pectral properties of three-dimensional d ynamical mod- 
els dFortnev et al.ll2006al:IShowman et alj|20 08). and now in pi- 
oneering the coupling of non-gray radiative transfer to three- 
dimensional dynamics (Showman et al. 2009). While these 3D 
models are being tested by observations of secondary eclipses 
and orbital phase variations with Spitzer, it is also possible to 
test their predictions for chemical mixing ratios and upper at- 
mospheric temperatures via transmission spectra observations. 
The main focus of this paper is to describe the results of the cal- 
culation of the transmission spectra of these three-dimensional 
models. 

The paper is organized as follows. In §2 we describe our 
methods, while in §3 we show results for the transmission spec- 
trum of simple isothermal models to investigate the role of sur- 
face gravity and enhanced metallicity. In §4 we describe our 
three-dimensional atmosphere models, while in §5 we review 
the transmission spectrum data for planet HD 189733b, and 
present our models, while §6 is the corresponding section for 
planet HD 209458b. §7 focuses on the observability of chem- 
ical abundance gradients in hot Jupiter atmospheres while §8 
gives our conclusions. 

2. METHODS: TRANSMISSION SPECTRUM 

We have previously developed a code to compute the 
transmission spectrum of planetary atmospheres. The first 
generation of the code, which used one-dimensional atmo- 
spheric pressure-te mperature (P-T) profiles, is described in 
Hub bard et al.l d2001|). In addition to ab sorption by atomic and 
molecular species. iHubbard et alj d200ll) also explored the ef- 
fects of refraction and a glow of photons around the planet's 
limb due to Rayleigh scattering. These additional effects were 
foun d to be n egligible for hot Jupiters, so we ignore them 
here. IFortnev et al. d2003l) investigated simple two-dimensional 



models of the atmosphere of HD 209458b, which included 
one planet-wide P-T profile, but 2D changes in the absorption 
cross-section of neutral Na, due to photoionization of Na as 
function of depth and angle from the substellar point. 

One can imagine a straight path through the planet's atmo- 
sphere, parallel to the star-planet-observer axis, at an impact 
parameter r from this axis. The gaseous optical depth tg, start- 
ing above the terminator and moving outward in one direction 
along this path, can be calculated via the equation: 

f°° r'dr'a(r')n(r') 

TG = l (r*-r^ ' 

where r' is the distance between the local location in the at- 
mosphere and the planetary center, n is the local number den- 
sity of molecules in the atmosphere, and a is the wavelength- 
dependent cross-section per molecule. We assume hydrostatic 
equilibrium with a gravitational acceleration that falls off with 
the inverse of the distance squaredQ In practice, the model at- 
mosphere is broken up into 1000 radius intervals, between the 
arbitrary base radius r and the maximum radius r + Ar, where 
Ar is 10-15xl0 3 km. The base radius, which is either at 10 
or 100 bars, is adjusted to yield the best fit to the observa- 
tions. Along the light ray path through the atmosphere parallel 
to the star-planet-observer axis, the local atmospheric density 
and opacity are also sampled at 1000 points. For absorption, the 
wavelength dependent cross-section is calculated based on the 
local chemical equilibrium mixing ratios, while the Rayleigh 
scattering cross-section is 

87r 3 (2+i/)V 

CTR= 3A 4 n 2 ' (2) 
where v is the refractivity (refractive index minus 1) of the gas. 
Since v <C 1 and v oc «, or is a function only of A and the gas 
composition. 

Here we define the wavelength-dependent transit radius as 
the rad ius where the total slant opt i cal dep th reaches 0.56, fol- 
lowing lLecavelier des Etangs et al.l ( 1200 8b). We match the ana- 
lytical results of these authors for the transit radius as a function 
of wavelength for a constant gravity, pure Rayleigh scattering 
atmosphere. 

Now we have modified and generalized the code to compute 
transmission spectra for a three-dimensional temperature struc- 
ture on an arbitrary latitude and longitude grid. We define a grid 
based on the star-planet-Earth axis, rather than the planet's rota- 
tion axis. The planet's substellar point is where the stellar zenith 
angle is zero. The zenith angle increases to 90 at the terminator 
and 180 and the antistellar point, which faces Earth during the 
transit. We additionally define a azimuth whose poles are the 
substellar and antistellar points. See Figure Q] The true north 
pole of the planet (defined by the rotation axis) runs through 
the azimuth=0 (prime) meridian. Azimuth values increase to- 
wards the trailing hemisphere. Using this new grid, we can first 
simply calculate the wavelength-dependent transmission of the 
planet at an azimuth, as a function of zenith angle. In prac- 
tice only zenith angles from ~ 70 to 110 degrees are sampled 
on the limb. Once the planet radius as a function of azimuth 
is determined, we then average these radii around the azimuths 
of the planet, to find the mean planetary wavelength-dependent 
radius. 

For any particular column of atmosphere, which is at a par- 
ticular azimuth/zenith angle location, hydrostatic equilibrium is 

5 While the transmission spectrum calculation allows for a decreasing surface gravity with height, the calculation the initial P-T profiles does not. The change in 
gravity with height is a minor effect, so this does not lead to a significant inconsistency. 
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assumed, and we use the given local P-T profile to interpolate 
in a pre-tabulated chemical equilibrium and opacity grid that 
exte nds out to 1 /ibar. The equilib r ium chemistry mixing ra- 
tios dLodders & Feglevll2002l l2006t lLodders||2009t) are paired 
with the opacity database dFreedman et al.l [20081) to generate 
pressure-, temperature-, and wavelength-dependent absorption 
cross-sections that are used for that particular column. For a 
different column of atmosphere, with a different P-T profile, 
local pressures and temperatures will yield different mixing ra- 
tios and wavelength-dependent cross-sections. 

3. CALCULATIONS FOR SIMPLE ATMOSPHERES 

Our survey of the literature finds that there has been no 
systematic exploration of the change of transmission spec- 
tra with planetary T e g for close-in planets. Other authors 
have previously ex plored reflection and e mission spectra (e.g . 



Marl ev et all 1 1 999b ; ISudarskv et all 120001: iBarmanet all 120011: 



Suda rskv et al . 2003; Marle v et all 120071) . but here we briefly 



examine transmission spectra for these atmospheres. To main- 
tain simplicity, we have investigated isothermal model atmo- 
spheres at 2500, 2000, 1500, 1000, and 500 K at two surface 
gravities, g = 10 and 50 m s" 2 . The chemistry and opacity 
databases are the same as those used for the 3D model atmo- 
spheres. The results of these calculations are shown in Figures 
[2] and [3] for the low and high gravity cases, respectively. All 
model planets have a base radius of 1.25 Rs at a pressure of 10 
bars. 

The striking impact of planetary surface gravity is clear from 
a comparison between these plots. While the absorption fea- 
tures (where the radius is larger) are the same, the dynamic 
range is far larger in the lower gravity planets. As previo usly 
discussed by other authors (e.g. lSeager & Sasselovl uOOO). the 
region of the planetary atmosphere probed in transmission is 
an annulus with an area proportional to the atmospheric scale 
height. Since lower gravity planets have larger scale heights, 
they are much better targets for transmission spectroscopy. 

In these calculations we assume chemical equilibrium. In 
the optical, this yields absorption bands of TiO and VO that are 
quite strong at 2000 and 2500 K, and weak at 1500 K. Sodium 
and potassium doublets are clearly strong at 1000 and 1500 K. 
In the 500 K model only bands of CH4 and H2O remain. In the 
blue-optical, Rayleigh scattering can be seen as a monotonic in- 
crease in radius towards the blue, although in the hottest models 
TiO absorption swamps most of this effect until blue-ward of 
400 nm. The slope of this Rayleigh scattering region is indica- 
tive of the temperature, and beco mes steeper (a more negative 
slope) as temperature increases dLecavelier des Etangs et all 
2008a). Also importantly the TiO/VO cross-sections fall off 
dramatically blue-ward of 0.4 /im. The detection of a dramatic 
change in radius here would likely indicate TiO absorption. In 
the infrared, bands of H2O and CO are strong in the models 
from 1500-2500 K, while CH 4 replaces CO in the 500 and 
1000 K models. Figures [2] and [3] should serve as more accu- 
rate replacements for the tr ansmission spectra of pM and pL 
class planets in Figure 1 1 of iFortnev et al.l ((2008). 

We next turn our attention to the effects of atmospheric 
metallicity for H/He dominated atmosphere s, as expected for 
close- in Jupiter- and Neptune-class planets. iMiller-Ricci et al.l 
(2009) have previously investigated the transmission spectra 
of 1-10 M ffi planets with hydrogen-rich and hydrogen-poor at- 
mospheres, and found that very large mean molecular weight 
(MMW) of hydrogen-poor atmospheres makes them extremely 



difficult targets for transmission spectroscopy, due to the much 
reduced scale height. In Figure|4]we investigate hydrogen-rich 
1000 K atmospheres from 1 x to 30 x solar metallicity. We find 
that while the MMW climbs from 2.35 to 2.97 in this case, the 
transmission spectra signatures are, to first order, not hampered. 

Figure|4] shows a lx model in black and 30 x model, with 
this same MMW, in thin green. Radius values at all wave- 
lengths increase due to the higher opacity everywhere. The 
thick green curve, with the correct MMW, brings these radii 
down to values similar to the lx model. Of interest at 1000 
K is that the atmosphere becomes much more CO and CO2 
dominated as the metallicity increases, since CO becomes more 
abundant linearly, and CO2 quadratically, w ith metallicity (e.g. 
lLodders & Feglevll2002t IZahnle et al.ll2009bl) . The strong C0 2 
bands at 4.4 and 15 /im are excellent and clear probes of at- 
mospheric metallicity. In the future, achieving spectra between 
4 and 5 /im will particularly useful, as CO is also very strong 
from 4.5 to 5 /im. The IRAC 4.5 /im band is so wide that it 
encompasses the CO and CO2 absorption features. 

4. METHODS: ATMOSPHERE MODELS 

We have run three-dimensional numerical simulations using 
the Substellar and Plane tary Atmospheri c Radiation and Circu- 
lation (SPAR C) model dShowman et alj|2009l) . which couples 
the MITgcm dAdcroft et al.ll2004l) to the plane-parallel, m ulti- 
stream radiative-transfer code o flMarlev & McK ay (1999). The 
MITgcm is a state-of-the-art atmospheric and oceanic circula- 
tion model that solves the global primitive equations, which are 
the standard equations for the large-scale dynamics in planetary 
atmospheres, in spherical geometry using pressure as a vertical 
coordinate. The primitive equations are generally appropriate 
when the atmosphere is stably stratified and when the horizon- 
tal/vertical aspect ratio of the atmospheric circulation is large, 
as is expected to be true on many hot Jupiters. The equations are 
discretized using finite-volume methods, which cast the equa- 
tions into a flux form that has good conservation properties, and 
adopt a "cubed sphere" grid rather than the more commonly 
used longitude-latitude grid, which allows longer time-steps 
and increases the accuracy near the poles. After completion of 
a simulation, the output can be interpolated ont o a regular lat- 
itude/longitude grid. See lShowman et ail d2009l) for additional 
details. 

Many previous atmospheric circulation models have forced 
the dynamics using highly idealized prescriptions for the ther- 
modynamic h eating/cooling rate, such as a Newtonian relax- 
ation scheme dShowman & G uillot 200]3 ICooper & Showman] 
2005, I2006h IShowman et alj 120081: iM enou & Rauscheri 
200 9; iRauscher & Menoul 120091) or flux-limited diffusion 
dPobbs-Dixon & Linl 120081) . In contrast, the SPARC/MITgcm 
model instead determines the thermodynamic heating/cooling 
rate using the state-of-the- a rt no n-gray radiative transfer 
scheme of iMarlev & McKa v (1999). The three-dimensional 
temperature structure can be viewed as numerous one- 
dimensional vertical columns (one per grid-box on the hori- 
zontal grid). At each model time-step, the model calculates 
the upward and downward intensities versus wavelength for 
each such vertical column. To calculate the thermodynamic 
heating rate over the full 3D grid, the model then sums the 
intensities over wavelength and takes the negative of the ver- 
tical divergence of the resulting net flux; this heating/cooling 
rate (which generally gives heating on the dayside and cooling 
on the nightside) is then used in the thermodynamic energy 
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equation to drive the dynamics. 

The planet's thermal radiative transfer is calculated for wave- 
lengths from 0.26 to 300 /im, with incident stellar flux treated 
from 0.26 to 6 /xm, using a correlated-k method to treat the 
opacities dGoodv et alj|1989|). A full des c riptio n of the opac- 
ity database is given in iFreedman et a D (120081) . The three- 
dimensional simulations described here, as in Showman et al. 
(2009), use 30 wavelength bins (see Showman et al. 2009 for 
more details). For a given metallicity, the chemical mixing ra- 
tios and opacities are calculated assuming local chemical equi- 
librium as a function of local temperature and pressure, ac- 
counting properly for depletion of gas-phase s pecies via con- 
densation, follow K. Lodders and collaborators ( Loddersl 1 999 ; 
Lodders&Feglevl l2002t lLoddersI 120021: IVisscher et aiT T2006; 



Loddersll2 009). Note, however, that we include only gas-phase 



served secondary eclipse photometry, spectra, and lightcurves. 
The main goals of this work is to further constrain the 3D model 
atmospheres in a region of the atmosphere that is distinct from 
that probed at secondary eclipse, as well as to point to future 
observations that may help us better understand these planetary 
atmospheres. 

5. HD 189733b 

Since its detection by Bou chv et alj d2005l) . planet HD 
189733b has lept to the forefront of exoplanetary characteri- 
zation. It features a large planet-to-star radius ratio and bright 
parent star, which enables detailed study of its atmosphere via 
observations of both its emission and transmission spectrum. 
For t he transmission spectru m, optical da t a have been published 
bv lPont et all d2008h andlRedfield et all (120 08), near infrared- 
data by Swain et alJ (120081) and lSing et al.1 (120091). a nd 3.6 and 
5.8 urn IRAC m i d-infra red data by Beaul ieu et al.l (2008) and 
Ehre nreich et al.1 d2007l) — th is Spitzer data was then reduced 
again bv lDesertetal.ld2009l) . who also added new 4.5 and 8.0 
/im points. 

Soon after the discovery of HD 189733b, iFortnev et al.l 
(2006b) computed the first models of its atmosphere. They 
found that the planet was considerably cooler than that of HD 
209458b, and was likely more simila r to TrES- 1 , which had jus t 
been observe d in secondary eclipse dCharbonneau et ai] l2005). 
In particular, IFortnev et al.l d2006bl) noted the strong methane 
absorption features in their ID model, which was meant to rep- 
resent a planet-wi de average P-T prof ile. This prediction was 
confirmed with the IS wain et alj d2008) detection of methane in 

the planet's transmission spectrum. 

In Figure |5] we compare an update of the IFortnev et al.l 
(2006b) ID profile (in black) to average profiles (averaged over 
all azimuths) from the 3D model, as a function of zenith angle, 
near the terminatorQ The ID model is a planet-wide average 
profile that assumes that the incident flux is first cut in half (for 
a day-night average) and then that the incident flux arrives from 
fi=0.5 (meaning th e incident angle is 60° ), which is the day- 
side average (see Fortnev & Marlev 2007). As previously sug- 
gested, the ID profile is a good approximation to the true limb 
profiles from the 3D model. Although the CO/CH4 ratio is al- 
ways > 1 for these mean profiles, given chemical equilibrium 
and the location of the profiles in P-T-space, we should see 
CH4 absorption features, in addition to CO, in the transmission 
spectrum, since CH4 is abundant. The equilibrium chemistry 
mixing ratios, for C-N-O molecules, for the ID HD 189733b 
model are shown in Figure [6] 

The transmission spectrum of ID and 3D solar-metallicity 
chemical equilibrium models of HD 189733b are shown in 
Figure|7] The published observations are also shown. The right 
y-axis is cast as the pressure at the terminator, for the planetary 
radius that would cause the observed absorption depth. The first 
thing to note is the similarity between the two model spectra. 
Since the ID profile is a reasonable proxy for the 3D profiles, 
the transmission spectra look very similar. 

The second thing to notice is that the fit to som e of the data 
are not particularly good. The iPont et ail d2008l) optical data 
may well be due to a haze layer, perhaps due to small en- 
statite grains dLecavelier des Etangs et al ]|2008a). This is a rea- 
sonable suggestion, given that silicate grains condense around 
these temperatures. However, it is not clear why similar ob- 

6 We note that these azimuthally averaged profiles are not used in the actual transmission spectra calculations. We perform the transmission spectrum calculations on 
a column-by-column basis, as a function of position, in the atmosphere. The azimuthally averaged profiles merely show average conditions at the planet's limb. 



opacity; any opacity from clouds/hazes is neglected. 

Local chemical equilibrium is likely a simplification for 
these highly irradiat e d, dy namic atmospheres. For instance 
ICooper & Showmanl d2006l) have previously shown that car- 
bon chemistry may be out of equilibrium, with CO favored 
at the ex pense of CH4, as is known for brown dwarfs of sim- 
ilar T„ff dSaumon et al.ll200l 120061: [Hubenv & Burrows! 120071 
iGeballeet al.ll2009l) . This process is also well-studied in the 
solar system (e.g IPriiin & Barshavl Il977t lYunget all [l988; 
iBezard et alj [2002). Alternatively, chem ical mixing ratios 
may be modified b y photochemistry (e.g. iLiang et al] [2004; 
Zah nle et al.ll2009bl) . The importance of photochemistry at the 
particular location that we are interested in here, the termina- 
tor, is not clear. Given these caveats, for HD 189733b in par- 
ticular we will investigate the transmission spectra assuming 
both chemical equilibrium and a fixed CO/CH4 ratio, follow- 
ing |CoopeE&^Showman!| d2006l) . Given the importance of CO2 
at high metallicity, future work on nonequilibrium chemistry 
should include the effect of this molecule as well. 

Here, we analyze SPARC/MITgcm simulations of HD 
189733b and HD 209458b that adopt solar metallicity and use 
a horizontal resolution of C32, approximately equivalent to a 
standard global grid of 128 longitudes x 64 latitudes. (Here is 
interpolated onto 64 zenith angles and 128 azimuths.) The sim- 
ulations use 53 vertical layers, with a base at 200 bar and a top 
at 2 fibzr (see Showman et al. 2009). Our HD 189733b case ex- 
cludes TiO and VO from the opacity calculation, while the HD 
209458b case includes both these species due to their possible 
relevance in producing the stratosphere inferred on that planet 
dKnutson et al.l|2008l; iBurrows et al.]|2007[: IFortnev etal] |2008; 
IShowman et all l2009tlSpiegel et al.H2009l) . Our HD 209458b 
simulation is identical to that in Showman et al. (2009) (see 
their Figs. 16-21); our HD 189733b simulations is analogous 
to their nominal HD 189733b case (see Figs. 4-10 in Showman 
et al. 2009) except that it includes more vertical layers and ex- 
tends to 2 /ibar rather than 0.2 mbar. Both simulations assume 
tidal locking of the planets — i.e., a permanent day and night 
hemi sphere. 

In Showma n et al.l d2009l) we considered the thermal emis- 
sion of these planets as a function of orbital phase, post- 
processing the P-T grids to compute spectra as a function 
of or bital phase, using the methods outlined in IFortnev et al.l 
( 2006a). These spectra can also be integrated over Spitzer band- 
passes to compute thermal emission light curves. These model 
spectra and lightcurves can then be directly compared to ob- 
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scuring silicates wouldn't effect HD 209458b, although it is 
possible that its atmosphere is hot enough for enstatite parti- 
cles to vaporize in the upper atmosphere. (Compare Figure 5 
to Figure 8.) Abundant silicate particles at such a great alti- 
tude (pressures of ^0.5 to 50 /ibar to explain the observations) 
may indicate vigorous vertical mixing above a cloud base at 
^100+ bar. Alternatively, the obscuring material could poten- 
tially be somethin g other than s ilicates — perhaps a methane- 
derived soot (Zahnle et al. 20 09a) . 

In the near-infrared, we match well the location of absorp- 
tion features found by ISwain et"a D <H)08), 

but not their am- 
plitude. Both models are normalized to the local radius mini- 
mum at 2 /Ltm. What limits the amplitude of these near IR fea- 
tures, which are predominantly due to water vapor absorption, 
is a strong contribution fro m H2 collision induced absorption 
(CIA) (e.g. lBorvsowll2002h . The model can not reach a deeper 
minimum in absorption, especially at 2.1 /im, because of th e 
CIA opa city (see, e.g. Figu re 15 from lSharp & Burrowsl l2007). 
Recently Sing et alTd2009l) have analyzed transit photometry 
in two narrow bandpasses, over this same wavelength range, 
and in the b l uer ch annel do not reproduce the depth found by 
ISwain et"al] ((2008). In the mid-infrared, it is p otentially prob- 
lemati c to make detailed co mparisons, as the iBeaulieu et alj 
(2008) and lDesertet ail 12009) reductions of the same data sets 
do not agree. This highlights that transmission spectrum obser- 
vations are still very much on the edge of what can be done, 
so choices made in data reduction have direct consequences on 
the interpretation. In addition, for the HD 189733 system in par- 
ticular, starspots can in principle make observed eclipse depths 
change from month to month, which may effect the comparison 
of models to multiple datas ets taken at different times. 

If we try to fit only the iPesert et al.1 J2009) data set from 
3-10 /im, we basically concur with their findings that a large 
CO abundance could be required, given the large radius in the 
4.5 /jm band. Alternately, however, CO2 could be much more 
important than CO. While the mixing ratio of CO increases lin- 
early with metallicity the abundance of CO2 increases quadrati- 
cally (e.g. lZahnle et al.l2009bl) . Our fit to the IRAC 4.5 ^m data 
is much better at high metallicity, as a strong CO2 band at 4.4 
/im increases the planetary radius, along with CO, in this band. 
In Figure|8]we compare transmission spectra calculated at 1, 3, 
10, and 30 x solar, for the ID model, wi th the corresponding 
correct MMWs, to the IPesert et all (120091) data. We also over- 
plot the band-averaged model calculations, and find that we can 
fit near the la error bars with the 30 x solar model. Alterna- 
tively, it is likely possible to fit with only a solar abundance of 
H2O, but a large CO or CO 2 abundance. I t is n ot possible for 
for us to achieve a fit to the Bea ulieu et al.l (120081) mid-IR data 
points plotted in Figure [7] with our P-T profile, unless we artifi- 
cially decrease the planet's gravity by a factor of ~ 5, to expand 
the scale height by this same factor. 

These fits bring up a modeling issue that we have not been 
able to resolve. In general, we are unable to match the large 
variation in transit radius of the models of Tinetti and col- 
leagues. This is true for th eir optical spectrum of HD 189733b 
shown in lPont et al.1 d2008l) . as well as th e mid-infrared spectr a 
for HD 189733b and HP 209458b fromlTinetti et all d2007al) . 
even when we use the iTinetti et al.l (l2007al) P-T profiles. Of 
course even when using the same profiles, the chemical mix- 
ing ratios may be slightly different between our calculations 
and theirs, although that may not explain the dramatic dif- 
ferences. Pifferent opacity databases may be more impo- 



rant, as high opacity directly corresponds to large radii and 
low opacity to small radii. We have had bet t er success 
in comparisons with other auth ors dBarmanl 120071; ISing et al] 
20085 iMiiler-Ricci etaill2009l) . and we match exactly the an- 
alytic radius vs . wavelength relation for a Rayle igh scatter- 
ing atmosphere {Lecavelier des Etangs et al . 2008b). Recently, 
Mad husudhan & Seagerl d2009l) investigated a large range of 
compositions for HP 189733b, and the amplitude of their ab- 
sorption features may be intermediate between those of Tinnetti 
et al. and our own. Of course, data is the true constraint, and 
Tin etti and colleagues ha ve shown excel l ent ag reement with 
the lSwain et al.l (120081) and lBeaulieu et all 12009) observations. 
Good agreement between models with similar inputs is impor- 
tant for constraining model atmospheres and for deriving chem- 
ical mixing ratios in observed planets, so this issue should be 
addressed with some care in the near future. 

6. HP 209458b 

Significant attention has also been paid to HP 209458b, 
which for many years was the only hot Jupiter amenable 
to atmospheric cha racterization. Neut r al ato mic sodium was 
first discovered by Charbonneau et al] (|2000|). and later con- 
firmed by other authors dSnellen et alj|2008t ISing et af] |2008a: 



Langland-^hulaetalJ ll>09). There seems to be little agree- 
ment as to the abundance of Na, and it could be time- variable. 
Bar man et al .1 (120021) found that non-LTE level populations of 
Na could contribute t o a w eaker absorption feature at 589 
nm, and lFortnev et al.l ([2003 ) found that photoionization of Na 
c ould be important. 

ISing et al.l (l2008al) used the data of lKnutson et al.l (l2007bl) 
to produc e an optic a l transm ission spectrum from 300 to 770 
nm. The ISing et al] d20~0 8a) data show absorption via Na, po- 
tentially Rayleigh scattering (likely due to H2) in the blue, and 
perhaps additional abso rption, due to ano t her absorbe r , which 
may be TiO and VO dP6sert et al]|2008l) . ISing et al] (|2008b) 
suggest that the Na feature shape is best matched by an at- 
mosphere enhanced in Na at lower pressures and depleted at 
lower pressures, due either to c ondensation of Na int o solids, or 
photoionization. Very recently Beaulie u" et al] d2009) have pub- 
lished 4-band IRAC photometry of the transit of HP 209458b, 
which is the same mid-IR spec tr al coverage attained f or HP 
189733b bv lbesert et al] (120091) . IBeaulieu et all (120091) inter- 
pret the radius variation observations as b eing due essentially 
only to water vapor, similar to the work of Tinet ti et al] (|2007b) 
for HP 189733b. 

In Figure [9] we plot the azimuth-averaged P-T profiles for 
the 3P model of HP 209458b, near the planet's terminator re- 
gion, and com pare this to our IP planet-wide average model, 
updated from lFortnev et al] (120051) . It is clear that HP 209458b 
possesses a larger diversity in P-T profiles in this region than 
HP 189733b. Models on the day side of the terminator have 
a temperature inversion, due to heating of the stratosphere via 
absorption of optical flux by TiO. The night side of the termi- 
nator region lacks TiO. For the IP planet-wide average profile 
(black), we do not find a temperature inversion, although for a 
day-side average model (which uses a 2 x hi gher incident flux), 
we do find an inversion (iFortnev et al]|2008l) . 

The comparison of the IP and 3P models with the optical 
data o f ISing etalJd2008al) in Figure[10]dramatically shows the 
difference in transmission spectra that these temperatures in- 
duce. The IP model, which lacks any appreciable amount of 
TiO or VO, generally matches the Na absorption features as 
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well as the radius increase in the blue, here due to Rayleigh 
scattering off of H2. The thin core of the Na absorption fea- 
ture is too strong in our model, which bolsters the finding of 
Sing et alJ (1200 8a) that Na is depleted in the upper atmosphere. 
The 3D model also shows the broad, pressure-broadened ab- 
sorption of Na and K, with numerous weak TiO/VO bands at 
some wavelengths. These bands are due to a significant deple- 
tion of TiO/VO compared to solar abundances. As shown in 
Figures [2] and [3] a solar complement of TiO/VO leads to a dra- 
matically larger radius. The 3D model looks somewhat similar 
to the 1500 K (green) models in these earlier figures. 

This ba sic picture of deple ted TiO/VO is similar to the sug- 
gestion bv lDesert et al.l d2008l) of weak TiO/VO bands red-ward 
of the Na doublet, but we do not achieve as good of a fit, since 
we do not leave the molecular abundances as free parameters. 
A small c hange in temperatur e could well lead to a fit similar 
to that of iDesert etafl d2008l) . ISpiegel et ail d2009t) have re- 
cently highlighted the myriad difficulties in modeling the verti- 
cal distribution of TiO in a ID model. It seems quite plausible 
that TiO, should it exist on the day side of HD 209458b, and 
cause its temperature inversion, could be quite depleted at the 
planet's terminator, where temperatures are colder. However, 
the 3D model does not repro duce the observed mid IR sec- 
ondary eclipse photometry (see lShowman et aTl 2009) or trans- 
mission spectrum very well. 

To further explore the terminator of HD 209458b we also 
plot the transmission spectra contributions of the day and night 
hemispheres separately, in FigureQT] The night hemisphere 
(blue) shows no TiO/VO, while the day hemisphere (red) shows 
the weak TiO/VO bands. The planet-wide spectrum (black) 
shows that it is the day side contribution that dominates the 
transmission spectrum. This is for two reasons: First, the scale 
height is larger on the warmer day side. This "puffing up" 
masks the absorption features formed on the night side. This 
will likely happen for any planet. Second, for HD 209458b in 
the optical in particular, the opacity on the day side is far larger 
than that on the night side, which further masks the contribution 
from the night side. 

We next compare our models to the recent 4-band IR AC 
data on the transit of HD 209458b dBeaulieu et al.ll2009t) . In 
Figure[T2]we compare the 3D model at solar metallicity, along 
with the ID model at 1 x and 10 x solar, to the mid-IR data set. 
Here the base 10-bar radius is chosen to give a good match to 
the data in the optical for the ID models, and this same base ra- 
dius is used for the 3D model. While the solar metallicity mod- 
els generally reproduce the slope towards larger radii at longer 
wavelength, no ne of the models reprod uce the large radius vari- 
ations found by Beauli euet al.l (120091) . Together with the pre- 
viously published poor fit of the HD 209458b 3D model to the 
secondary eclipse data (Showma n et al.ll2009l) . this may indi- 
cate that the model does not reproduce the true mid-infrared 
opacity of the planet. We also note that while the ID and 
3D models have dramatically different optical transit radii, due 
to the temperature-sensitive TiO/VO abundances, the same is 
not true in the mid-infrared, where both the ID and 3D solar- 
metallicity models are dominated my H2O and CO opacity, and 
predict very similar radii. 

7. SPECTRAL DIFFERENCES OF THE TRAILING AND 
LEADING HEMISPHERES 

In the current era it is extremely difficult for observers to even 
achieve data on hot Jupiter atmospheres. However, in the future 



this will not always be the case. We are quite optimistic about 
the era of the James Webb Space Telescope, which is sched- 
uled to launch in 2014. It will then be possible to derive more 
detailed information about these planets' atmospheres. In par- 
ticular, transmission spectra of a planet's leading hemisphere 
(probed during the first half of transit ingress) and trailing hemi- 
sphere (during the last half of transit egress) may allow for 
detection of spectral differences on different locations of the 
planet. This could be complementary to eclipse mapping during 
the secondary eclipse ingre ss and egress phases (Will iams et al.l 
2006: lRauscher et alJl2007h . 

It is a general outcome of our published 3D hot Jupiter 
simulations that wide low-latitude west-to-east winds persist 
on arbitrarily long time scales. The first observational conse- 
quence is that the hottest point of the planet is blown down- 
stream from the substella r point (|Showman & Guillol 120021 : 
Coo per & Showman] 120051: I Showman et alj|2009t). w hich was 
confirmed for HD 189733b bv lKnutson et all d2007al) . The ef- 
fect is more dramatic in simulations of HD 189733b than for 
HD 209458b, o wing mostly to a longer radiative time con stant 
in HD 189733b dShowman et alJl200alR?rtnev et alj|200l . At 
the time of primary transit, this leads to the leading hemisphere 
being colder than the trailing hemisphere, due to the downwind 
displacement of both th e hottest and coldest points of the planet 
dShowman et alj|2009l) . 

This dichotomy could potentially be probed with transmis- 
sion spectra obtained during the transit ingress and egress 
phases. A transmission spectrum taken half-way into the time 
of ingress would sample only the leading hemisphere, and at 
egress, the trailing. The simple leading-trailing image is com- 
plicated by the transit impact parameter, b, however, as shown 
in Figure [13] As b increases from zero (central transit) to one 
(grazing), the edge of the parent star becomes gradually more 
inclined to the planet's rotation axis. ForHD 189733b this leads 
to 38.8° of trailing hemisphere (near the pole) on the leading 
edge, and for HD 209458b, 30.9°. In general this would serve 
to mute differences between the hemispheres. 

The differences in the spectra of the hemispheres can be seen 
most clearly in HD 189733b which is a relatively cool hot 
Jupiter. Its planet-wide mean P-T profile is close to the bound- 
ary wh ere CO and CH4 have equal abundances dFortnev et al.l 
2006b), given chemical equilibrium. We find that if this equi- 
librium holds, the night side s hould be CH4 domina ted, while 
the day side is CO dominated dShowman et al 1 120091 We have 
computed separate transmission spectra of the leading and trail- 
ing hemispheres (including the impact parameter effect) which 
are shown in Figure [14] The differences in temperature lead 
to dramatic differences in spectra. The leading (cooler) hemi- 
sphere has a smaller scale height, which leads to a smaller tran- 
sit radius at most wavelengths. However, in CH4 bands, absorp- 
tion is much stronger, which can lead to radii larger than that 
of the trailing hemisphere. The dramatic differences in spectra 
may be within reach of JWST, if multiple transits are observed 
(T. Greene, personal communication). 

These CO/CH4 absorption features could also be a direct 
probe of non-equilibrium chemistry in the atmosphere of HD 
189733b, since vigorous vertical mixing would tend to ho- 
mogenize t he CO/CH4 ratio everywhe re in the planet's atmo- 
sphere (see Cooper & Showman 2006, and references therein). 
In order to simulate how a fixed CO/CH4 ratio, independent of 
height and location, would effect these spectra, we have calcu- 
lated transmission spectra with mixing ratios given by the "cold 
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model" of lCooper & Showmanl (12006). which is a good approx- 
imation for the temperatures of HD 1 89733b. That model yields 
a CO/CH4 ratio of 4.5/1, meaning that CO dominates, but both 
molecules are abundant and should show strong absorption fea- 
tures. The results are shown in Figure[15] The spectra of the 
hemispheres are essentially identical, except for the larger scale 
height on the trailing hemisphere. Absorption features of both 
CO and CH4 can be seen. 

For HD 209458b, which is significantly hotter than HD 
189733b, we find smaller differences between the leading and 
trailing hemispheres. The warmer temperatures lead both hemi- 
spheres to be CO-dominated, with CH4 absorption features ab- 
sent. In addition, the shorter radiative time constant in the pM 
class atmos phere leads to smalle r shift in the hot spot down- 
stream (see lShowman et al] |2b09). Figure [16] shows difference 
in the hemispheres are mostly due to small scale height differ- 
ences, with some subtle effects on TiO absorption bands, due 
to its more deplete d abundance on the cooler (leading) hemi- 
sphere. Previously llro et aL (2005) showed a similar effect for 
a 3D model of HD 209458b that included radiative transfer and 
parameterized winds. This model lacked TiO, and they found 
a strong Na absorption feature at 589 nm on the hotter trail- 
ing hemisphere, but a weaker Na feature on the cooler leading 
hemisphere, due to the condensation of Na into clouds. They 
did not investigate transmission spectra at other wavelengths. 

8. CONCLUSIONS 

We have investigated several aspects of the transmission 
spectra of transiting planets. As previously known, high tem- 
perature and low surface gravity planets yield the best prospects 
for obtaining transmission spectra. We have investigated the 
spectra for a model planet at an atmospheric metallicity of up 
to 30 x solar, which is likely representative of Neptune-class 
exoplanets, and demonstrated that enhanced MMW weight in 
these atmospheres does not significantly diminish the magni- 
tude of absorption features in these atmospheres. The small 
cross-sectional area of the atmospheres of these smaller planets 
is a greater challenge. 

We have calculated the transmission spectrum of 3D models 
of the atmospheres of HD 189733b and HD 209458b, and com- 
pared these results to those of the ID calculations. Differences 
between the 3D and ID models can be dramatic, or negligible, 
depending on the wavelength of observation, and the tempera- 
tures probed at the planet's terminator region. For the planet- 
wide 3D transmission spectrum of HD 189733b, we find it is 
well reproduced by the ID model. Our mode l constrains the 
postulated upper haze layer dPont et alj|2008l) to be optically 
thick at optica l wavelengths at microbar pressures. Following 
iFortnevI (120051) . we find that the slant optical depth at 1 /ibar 
of an absorbing layer that shares the gaseous scale height is 60 
times larger than the normal optical depth. It remains possi- 
ble that such an absorber could remain undetected in scattered 
light. 

We are able to reproduce the IRAC transmission spectrum 
points of Des ert et alj d2009l) with a large abundance of CO, 
and especially CO2, in these models due to a large supersolar 
metallicity. This is at odds with the determination of subsolar 
metallicity for this planet by comparison to near-i nfrared data 
dSwain et all 120081; iMadhusudhan & Seageril2009l). Th e pub- 



lished 3D model of HD 189733b from lShowman et alj ((2009) 
has been able to reproduce most of the secondary eclipse pho- 



tometry and spectrosco py, the broad features of the 8 and 24 fim 
half-orbit light curves dKnutson et aHl2007al |2009), and some 
features of the transmission spectrum (although an unmodeled 
haze apparently dominates in the optical and possibly the near 
IR). This gives us confidence that these models, which are still 
in their infancy, are on the right path towards being a realistic 
simulation of dynamic hot Jupiter atmospheres. 

For planet HD 209458b, the differences between 3D and ID 
models of the transmission spectrum are significant at optical 
wavelengths, due to the wide range of temperature probed in 
the upper atmosphere of HD 209458b at the terminator region. 
In the 3D model we find TiO in abundance on the day side 
dShowman et al.l2009l) and a reduced abundance on the planet's 
limb. This TiO masks absorp tion due to Na and K (which could 
well be photoionized — e.g. lFortnev et al.ll2003l). and we find a 
poor match to the transmission spectrum of ISing et alj d2008a). 
The higher temperatures and higher opacity of the day-side of 
the limb mask absorption features from the night side of the 
limb. The ID model, a planet-wide average, shows no TiO ab- 
sorption features, and is a better match to the optical data. The 
poor match to the transmissions spectru m, and to secondary 
eclipse photometry (Sh owman et alj2009f) indicates that the 3D 
model presented here and in Showman et al. (2009), which as- 
sumes chemical equilibrium and produces a day-side tempera- 
ture inversion using TiO and VO, is not a good representation 
of the atmosphere of HD 209458b. It is possible that the ab- 
sorbers causing the temperature inversion in this atmosphere 
do not have the same optical properties and/or vertical profile 
of the TiO gas employed in our model. Future 3D models will 
explore the effects of disequi librium compounds, such as sulfur 
species (Zahnle et al. 2009b), in producing thermal inversions. 

Transmission spectra will continue to be an important tool in 
characterizing the atmospheres of transiting planets, with HST 
in the UV to near IR, as well as with JWST in the red optical 
to mid IR. Probing the atmospheres with this method will con- 
tinue to be complementary to secondary eclipse spectra, since 
different regions of the atmosphere (in latitude, longitude, and 
height) are probed with each method. In addition, transmission 
spectra are less sensitive to the temperature gradi ent, while this 
is a m ajor complication for emission spectra (e.g. lFortnev et alj 
I2006al) . 

One could imagine more challenging levels of atmospheric 
characterization, such as the detection of wind speeds via 
the Dop pler sh ift of absorption lines in the planetary atmo- 
sphere dBrownl 1200 lb . This could potentially be done dur- 
ing the ingre s s and egress phases. However, as discussed by 
Spiegel et al. (2007), the characterization of transiting planet 
atmospheres at this level requires the consideration of a num- 
ber of additional effects, such as the planet's rotation rate and 
orbital motion, amongst others. We will consider this in the 
future. 
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FIG . 1 . — Diagram of our zenith angle and azimuth grid, that aids in the computation of 3D transmission spectra. The antistellar point faces the observer, while the 
substellar point face the star. Azimuths are in red, while zenith angles are in green. The terminator is shown in blue and the planet's true north pole is also labeled. 
The blue band and purple slice are meant to schematically show the region of the atmosphere that is being probed during the transit. 
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FIG. 2. — Planet radius vs. wavelength for isothermal model atmospheres with a gravitational acceleration of 10 m s 2 . Models are T=2500, 2000, 1500, 1000, 
and 500 K, from top to bottom and assume chemical equilibrium. Prominent absorption features, which appear as increases in the planetary radii, are labeled. 
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FIG. 3. — Planet radius vs. wavelength for isothermal model atmospheres with a gravitational acceleration of 50 m s 2 . Models are T=2500, 2000, 1500, 1000, 
and 500 K, from top to bottom and assume chemical equilibrium. 
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FIG. 4. — Planet radius vs. wavelength for the 1000 K isothermal model atmosphere at 4 values of metallicity. For the thick lines, the correct mean molecular 
weight is calculated for each metallicity. The thin green line is the 30 X solar model, but with an unphysically low mean molecular weight (that of the 1 X model) 
for comparison. 
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FIG. 5. — Atmospheric P-T profiles at the limb of HD 189733b, from the 3D model. Smaller zenith angles are towards the day side, and larger zenith ang les 
towards the night side. From day-side to night-side, these are red, orange, green, blue, and cyan. The ID planet-wide average P-T profile (Fortney et al. 2006b) is 
shown in black. The dashed black curve shows where CO and CH4 have equal abundances. CO is increasingly favored at higher temperatures. The condensation 
curve where TiO gas is lost to solid condensates is shown as a dotted black curve, as is the condensation curve for enstatite, MgSiOj . 
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FIG. 7. — HD 189733 in-transit absorption depth vs. wavelength for ID (black) and 3D (red) model a tmospheres, assuming chemical equilibrium. Computed 
absorp tion depths are (R B / R s) 2 - Many data sets are also shown. The optical (orange) is Pont et al. 1 2008) an d the near infrared (gTe en and purple) are Swa in et ail 
1 2008) and Sing et al. 12009), respectively. Further into the infrared, in dark blue at 3.6 and 5.8 /im are from Beaulieu et al. 12008), while the 8 /im point is from 
Knutson et al. (2007a), and the 24 fim point is from Knutson et al. 12009). The orange 8 /im point is from Agol et al. (2009). The four cyan IRAC points are from 
Desert et al. 12009). The data of Redfield et al. 12008) are of much higher spectral resolution that the other data sets, and are not plotted. The right y-axis shows the 
pressure at the terminator that corresponds to a given planet radius/absorption depth. 
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FIG. 8. — In-transit absorption depth vs. wavelength for I D models of HP 18 9733b at four metallicities, from 1 X to 30x solar. Computed absorption depths are 
(Rp/ Rs) 2 - We compare in detail to the mid-infrared data of Desert et al. ( 2009), from 3-10 [im. Band-averaged model absorption depths are shown as squares. A 
higher metallicity favors CO2 and CO, at the expense of CH4, yielding a better fit to the 4.5/im band in particular. 
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FIG. 9. — Atmospheric P-T profiles at the limb of HD 209458b, from the 3D model. Smaller zenith angles are towards the day side, a nd larger zenith ang les 
towards the night side. From day-side to night-side, these are red, orange, green, blue, and cyan. The ID planet-wide average P-T profile iFortney et al. 2005) is 
shown in black. The dashed black curve shows where CO and CH4 have equal abundances. CO is increasingly favored at higher temperatures. The condensation 
curve where TiO gas is lost to solid condensates is shown as a dotted black curve, as is the condensation curve for enstatite, MgSiO^ . There is clearly a larger spread 
in the temperatures at the limb of HD 209458b than for HD 1 89733b. 
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FIG. 10.— Transit radii of HD 209458b ID (black) and 3D (red) models, compared to the data of ISing et alj (2008a), in green. The models, which assume 
chemical equilibrium, have been binned onto the Sing et al. wavelength grid. The ID model shows Rayleigh scattering in the blue, along with a broad Na absorption 
feature at 589 nm, that is stronger than observed. The 3D model shows wea k bands due to TiO and VO, that are still much stronger than observed. The data are 
plotted assuming a parent star radius of 1 . 125 Kg, after Knutson et al. 1 20 07bT) . 
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FIG. 1 1 . — Transmission spectrum contributions of the day (red) and night (blue) hemispheres of the planet. The full-planet model is shown in black. The larger 
scale height and higher optical opacity of the day side lead it to dominate the planet-wide transmission spectrum. 
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FIG. 12.— Transit ra dius of HD 209458b ID and 3D models, compared to the IRAC data of iBeaulieu et all <2009l) and the mean 24 /im MIPS data point of 
Richardson et al. (2006) and H. Knutson (pers. comm.). The 3D model (which is 1 X solar) is shown in purple, while ID models at 1 X and lOx solar metallicity are 
shown in red and orange respectively. The base radius is chosen to give a good fit to the optical observati ons. Ba nd-averaged model radii over the IRAC bandpasses 
are shown in squares. The data are plotted assuming a parent star radius of 1.125 Rs, after Knutson et al. (2007b). Water vapor and CO features dominate the mid-IR 
spectrum, as in HD 189733b, while the 4.4 fim feature in the 10 X model is due to CCV 
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FIG. 13. — Schematic view of ingress at a non-zero impact parameter. Half-way through ingress, one samples mostly the leading hemisphere, but also some of 
the trailing hemisphere, near the pole. The planet's rotation axis is the long-dashed line, while the planet's orbital path is the dotted line, with arrows. 
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FIG. 14. — Planet radius vs. wavelength for the leading and trailing hemispheres of a 3D simulation of planet HD 189733b, assuming chemical equilibrium. The 
leading hemisphere (blue) is cooler, while the trailing hemisphere (red) is warmer. Cooler regions have a smaller CO/CH4 ratio, which leads to changes in spectral 
features on the two hemispheres. The leading hemispheres shows much stronger CH4 absorption. 




FIG. 15. — Planet radius vs. wavelength for the leading and trailing hemispheres of a 3D simulation of planet HD 189733b, with a fixed CO/CH4 mixing ratio 
of 4.5/1, as could be expected from homogenization due to vertical mixing. The leading hemisphere (blue) is cooler, while the trailing hemisphere (red) is warmer. 
Differences between the hemispheres are due only to scale height differences, and the small temperature-dependent changes in opacity. 
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FIG. 16. — Planet radius vs. wavelength for of the leading and trailing hemispheres of a 3D simulation of planet HD 209458b, assuming chemical equilibrium. 
The leading hemisphere (blue) is cooler, while the trailing hemisphere (red) is warmer. Most the differences are due to the change in the scale height, while the 
TiO/VO band are stronger on the warmer, trailing, hemisphere. 



